Algebraic spin liquid as the mother of many competing orders 
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We study the properties of a class of two-dimensional interacting critical states - dubbed algebraic 
spin liquids - that can arise in two-dimensional quantum magnets. A particular example that we 
focus on is the staggered flux spin liquid, which plays a key role in some theories of underdoped 
cuprate superconductors. We show that the low-energy theory of such states has much higher sym- 
metry than the underlying microscopic spin system. This symmetry has remarkable consequences, 
leading in particular to the unification of a number of seemingly unrelated competing orders. The 
correlations of these orders - including, in the staggered flux state, the Neel vector and the order 
parameter for the columnar and box valence-bond solid states - all exhibit the same slow power-law 
decay. Implications for experiments in the pseudogap regime of the cuprates and for numerical 
calculations on model systems are discussed. 



I. INTRODUCTION 

In the effort to explain the still-mounting puzzles in 
many strongly correlated materials, one frequently in- 
voked idea is that of competing orders. Specifically, it 
is often appealing to contemplate the presence of rather 
slowly-varying fluctuations in two or more different or- 
der parameter degrees of freedom. In some cases these 
orders are not obviously related to one another - one oft- 
discussed example is antiferromagnetism and supercon- 
ductivity in the cuprate high-T c superconductors^ This 
kind of situation raises an important question: are the 
competing orders controlled, all together, by the univer- 
sal physics of a single phase or critical point? We can also 
turn this question on its head - rather than phcnomcno- 
logically introducing some set of slowly-fluctuating or- 
ders, we can take a somewhat more microscopic approach 
and look for possible quantum states of a given system. 
Then we can ask whether competing orders arise natu- 
rally in some such state. 

In this paper we shall follow this strategy and show 
that, somewhat surprisingly, this physics obtains within 
a certain spin liquid stated of two-dimensional electronic 
Mott insulators that has been suggested to play a key role 
in the underdoped cuprate superconductorsi^^*SiL& The 
particular spin liquid state we consider has been variously 
described as a "<i-wave" resonating valence bond (RVB) 
state or a staggered flux (sF) state. Here we will use 
the latter nomenclature and refer to it as the staggered 
flux state. It is important to note that the staggered flux 
spin liquid possesses no broken symmetries and is quite 
distinct from ordered states with a staggered pattern of 
orbital currents; instead, it is a specific incarnation of the 
RVB idea of Anderson^ 

Previous papers have shown that the sF spin liquid is 
an interacting critical state, and that it may be a sta- 
ble critical phasei^H^ - the spin correlations decay 
as a power of the distance with a universal exponent, 
and, while a description in terms of fractional 5=1/2 



A?' 




FIG. 1: Cartoon pictures of some of the slowly-varying com- 
peting orders within the staggered flux spin liquid state. 
These are the Neel state (A), and the columnar (B) and box 
(C) valence bond solids. The shaded regions denote those 
groups of spins that are most strongly combined into local 
singlets. Note that in the sF state these orders fluctuate in 
both space and time. These pictures describe the character of 
some of the important slowly- varying fluctuations, but should 
not be viewed as snapshots of the physics at the lattice scale, 
which may be quite complicated. 



spinons is natural, they do not behave as free quasiparti- 
cles even at asymptotically low energy. Furthermore, the 
dynamic critical exponent z = 1. Alternatively, the long- 
distance, low-energy properties are controlled by an in- 
teracting, conformally invariant fixed point. Such states 
were dubbed algebraic spin liquids (ASL) in Ref.l?]. Here 
we show that, remarkably, several competing orders are 
unified within the sF state by an emergent SU (4) symme- 
try, and all have the same slowly-varying long-distance 
correlations. 

Two of the competing orders are simply the the Neel 
vector, and the order parameter for the columnar and 
box valence-bond solid (VBS) states - cartoon pictures 
of these orders are shown in Fig. To be precise, con- 
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sider a model Mott insulator on the square lattice with 
one electron per site at zero temperature, and suppose 
the system has been tuned into the staggered flux spin 
liquid by increasing the magnetic frustration. This state 
should be present in the parameter space of this system 
at least as a multicritical point but potentially as a stable 
phase. A measurement of the Neel correlations will find 
the power-law decay 

(-l)(r*+r y ) {Sr . S(j) _ i/|r.|(i+n). (1) 

Remarkably the VBS correlations display the same be- 
havior; for example, 

(-l) r *{(S r+x ■ S r )(S x ■ S )) ~ l/|r|( 1+ "'. (2) 

This correlator measures the tendency of the system to 
order into the columnar dimer pattern shown in Fig. ^3. 
Furthermore, several other (more unusual) competing or- 
ders with the same power-law decay are present. These 
are the plaquette-centered spin at q = (it, it), the density 
of Skyrmions in the Neel vector at q = (jr, it), and a kind 
of triplet valence bond order that breaks spin rotations 
but not time reversal. Also, the uniform spin chirality 
exhibits slow power-law decay with an exponent that is 
likely the same as for the other orders. It is important 
to note that the decay may be quite slow and thus is po- 
tentially observable in experiments and numerical simu- 
lations - variational wavefunction studiesi&i^ of the Neel 
correlations provide the rough guess r\ w 0.5. 

This physics may have important consequences for the 
pseudogap regime of the underdoped cuprates. It has 
been suggested that this part of the phase diagram may 
be viewed as a doped sF spin liquidi 5 ' 6 ' 7 ' 8 In particular, 
in the spin sector, the rather high-temperature physics 
of the pseudogap should be essentially unchanged from 
the undoped sF state. The presence of slowly-fluctuating 
competing orders related by a large SU (4) symmetry 
opens up a new possibility for tests of this hypothesis. 
The simplest of these to probe is almost certainly the 
Neel vector, the fluctuations of which are directly mea- 
sured by magnetic neutron scattering at q = (it, it). Fur- 
thermore, the sF state is described by a critical theory, 
so if it is present this magnetic scattering should exhibit 
critical scaling. If this is found it will be important to 
think about whether the other competing orders related 
to the Neel vector by SU(A) symmetry can be directly 
probed. These points are discussed in more detail in 
Sec. I VII - readers not interested in following the more 
technical aspects of our results may wish to skip to this 
point. 

The sF state is just one member of a class of ASLs that 
also give rise to a variety of competing orders unified by 
a large emergent symmetry. Another state of particu- 
lar interest is the 7r-flux (7rF) spin liquid^, of an SU (4) 
Heisenberg-like model on the square lattice. Assaad has 
recently carried out quantum Monte Carlo simulations 
of this model; 1 — the results suggest that the irF state 
may have been observed, but further tests are desirable. 



The results of this paper can be tested numerically and 
should significantly aid the resolution of this issue. We 
give concrete suggestions along these lines in Sec. [V] 

From a formal point of view, the sF and ttF states can 
both be described at low energies by a field theory of 
fcrmionic spinons with massless Dirac dispersion, mini- 
mally coupled to a noncompact U(l) gauge field. This 
theory is often referred to as noncompact QED3. There 
is good evidence that it can flow to a conformally invari- 
ant, interacting fixed point, over which one has control in 
the limit of a large number (Nf) of fermion flavors^&i 7 . 
This fixed point is the algebraic spin liquid. Here Nf is 
the number of two-component Dirac fermion fields, and 
these can be rotated into one another by an SU(Nf) fla- 
vor symmetry - we have Nf = 4 for the sF state and 
Nf = 8 for the ttF state of the SU(A) spin model. It 
has been shown that, in the large-iV/ limit, all relevant 
perturbations to the ttF state are forbidden by symmetry 
and it is thus a stable phased The same conclusion is 
easily seen to hold for the sF state (see Sec. 1111 B|) . It 
is not known whether stability continues to hold for the 
interesting values of Nf, although the results of Ref. ITU 
suggest that the 7rF state is stable for Nf = 8. Even 
if the sF state is unstable, it should appear as a zero- 
temperature critical or multicritical point and may still 
be interesting. 

In the field theory, the competing orders discussed 
above arise as follows: In the simplest scenario, which 
is suggested by the 1/iVy expansion, the dominant corre- 
lations are those of an SU(Nf) adjoint N a and a scalar 
M - these are bilinears of the fermions. It is a simple 
matter to work out how these operators transform under 
the symmetries of the spin model, and to find symmetry- 
equivalent physical observables with the same transfor- 
mation properties; these quantities will all exhibit power- 
law correlations decaying as l/|r| 2Alv or l/|r| 2Aju , where 
A at, Am < 2 are the scaling dimensions of N a and M, 
respectively. (In fact An = Am to all orders in 1/Nf, 
although it is not clear whether this holds at finite Nf.) 
Note that 2A^ = 1 + rj. Both the Neel vector and the 
order parameter for the columnar and box VBS states 
are symmetry-equivalent to particular components of N a ; 
this gives rise to the power-law decay of Eqs. Q and 

We note that the structure of competing orders aris- 
ing from QED 3 has been discussed previously, from a 
rather different point of view, in a different physical 
context .IgiiiLSli Also, it was recently observed that Neel 
and VBS orders can be unified (at the mean-field level) 
by a chiral rotation at the irF saddle point i2i 

We now outline the rest of the paper. We review the 
description of the sF spin liquid in Sec. [H] Section Til Al 
discusses the route from the slave-fermion description of 
the Heisenberg model to the field theory, and Sec. Ill Bl 
reviews the use of the large- Nf expansion to control the 
sF fixed point. In Sec. lIII Al we discuss in detail the sym- 
metries of the sF state and their associated conserved 
currents. Section Till Bl extends the argument of Ref. IT2I 
for the stability of the ttF state at large Nf to the sF state 
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- the only significant difference is the presence of velocity 
anisotropy, which is dealt with in Refs. |52 and and 
Appendix|0 Our main result for the sF state is the iden- 
tification of the slowly-varying competing orders - this is 
discussed in Sec. IIV Al The same is done for some com- 
ponents of the conserved currents in Sec. IIV Bl In Sec. El 
we shift gears to discuss an analogous identification of 
competing orders for the irF state of an SU(4) Heisen- 
berg model. Finally, in Sec. IVII we discuss the prospects 
for observation of this physics in the cuprates, and we 
conclude in Sec. IVHI with a discussion of some of the 
issues raised by our results. 

II. DESCRIBING THE ALGEBRAIC SPIN 
LIQUID 

A. From the lattice to the continuum 

We begin by reviewing the description of the algebraic 
spin liquid fixed point, and the staggered flux state in 
particulari 5 ! 6 ! 7 ! 10 ! 11 ! 12 ! 23 ! 24 The starting point is the slave- 
fermion mean-field theory of the S = 1/2 Heisenberg 
model on the square lattice: 

H = J Y S r ■ S r > + ■ ■ ■ ■ (3) 

{rr') 

Here J > (antiferromagnetic exchange), and the ellipsis 
represents perturbations consistent with the symmetries, 
such as further neighbor frustrating exchanges, ring ex- 
change terms, and so on. We require that the Hamilto- 
nian be invariant under SU(2) spin rotations, time rever- 
sal, and the full space group of the square lattice. For- 
mally we may rewrite the spin as a bilinear of fermionic 
"spinon" operators 

s r = ^ftoPapUp- ( 4 ) 

Here a = 1,2, and f^ rl (/J 2 ) creates a spin-up (spin- 
down) fermion. This is an exact rewriting when com- 
bined with the local constraint fl a f ra = 1. Exploit- 
ing the well-known SU (2) gauge redundancy in spinon 
variables the spin-spin interaction is decoupled with 
an SU{2) gauge field residing on the links of the lat- 
tice. The mean-field saddle points arc then described by 
quadratic spinon Hamiltonians; depending on the struc- 
ture of the saddle point, the important low-energy fluc- 
tuations enter via an SU(2), U(l) or Z 2 gauge field min- 
imally coupled to the spinonsiSI 

At the mean-field level, the sF state is described by 
the Hamiltonian 

W "f = -E E [(ft+(-l) (r *- r ' ) A)/t a / p , a + H.c] J 

r£A r'n.n.r 

(5) 

where the first sum is over sites in the A sublattice, and 
the second is over the nearest neighbors of r. This de- 
scribes spinons hopping in a background staggered flux of 



$ = ±4 arctan(£/A), where the sign alternates from one 
sublattice of square plaquettes to the other. The appar- 
ent breaking of translation symmetry is a gauge artifact; 
the spinons transform under certain lattice symmetries 
with an additional SU(2) gauge transformation. Physi- 
cal operators are gauge invariant, so their transformation 
properties are unaffected and the saddle point possesses 
the full symmetry of the microscopic model. This situ- 
ation is summarized by saying that the spinons obey a 
projective symmetry groupie (PSG). The action of the 
PSG on the spinons is specified in detail in Appendix [Al 
The low-energy fluctuations about Eq. JSJ are encapsu- 
lated by a compact U(l) gauge field minimally coupled 
to the spinons. The full lattice Hamiltonian takes the 
form 

H sF = h E e lr> ~ K E C0S ( CUrl °>) ( 6 ) 
{rr') □ 

-E E [(it+(-l)^- r ^A)fi a e- ia ^f r , a + K.c.]. 

r£A r'n.n.r 

Here e rr i and a rr i are lattice vector fields: e is the elec- 
tric field and takes integer eigenvalues, while a, the vector 
potential, is a 27r-periodic phase. On the same link of the 
lattice, e and a satisfy the canonical commutation rela- 
tion [a, e] = i. The second term of Eq. 10 is a sum over 
square lattice plaquettes, and (curia) is the discrete line 
integral of the vector potential taken counterclockwise 
around the given plaquette. The Hamiltonian must be 
supplemented by the gauge constraint 

(dive) r +/t a / ra = l, (7) 

where (dive) r is the lattice divergence of the electric 
field. This gauge theory reduces exactly to the nearest- 
neighbor Heisenberg model in the limit K = and 
h/t — > 00; in this limit e = and the gauge constraint be- 
comes frafra = !• We will be interested instead in a limit 
where the mean-field theory is manifestly a good starting 
point and is valid up to intermediate length scales, so we 
consider K ^> t 3> h. (In both cases A w t.) The result- 
ing ASL fixed point will control the low-energy physics 
for some spin Hamiltonians of the general form of Eq. J2J . 
If the fixed point is stable, no fine-tuning should be nec- 
essary to access this part of parameter space, but the 
precise microscopic requirements are unknown. 

In the limit of interest, the gauge fluctuations are 
strongly suppressed by the large Maxwell term; that is, 
fluctuations in (curl a) at the scale of the lattice are very 
small. We can therefore first write a continuum theory 
of the long- wavelength, low-energy free-fermion excita- 
tions of Eq. |jSJ|, and then include the gauge fluctua- 
tions. The technical details are outlined in Appendix [Al 
the resulting low-energy theory consists of four massless 
two-component Dirac fermions minimally coupled to a 
noncompact U(l) gauge field. Microscopically, the gauge 
field is compact, which means physically that instanton 
configurations (magnetic monopoles) are allowed in the 
action. Therefore we should view the noncompact theory 
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as a point in the parameter space of the compact theory 
where all monopole fugacities have been tuned to zero 
- we shall be interested in expanding about this point. 
As discussed in Ref. El monopoles can (and must) be 
incorporated as perturbations. This is greatly aided by 
the observation that the absence of monopoles is pre- 
cisely equivalent to the presence of an emergent global 
£^(l)flux symmetry corresponding to the conservation of 
gauge flux, which is only conserved modulo 2tt in the 
compact theoryiSSiSS 

For simplicity of notation, it is convenient to suppress 
all fermion indices and work with the eight-component 
object 'J. We express matrices acting on ^ as tensor 
products of the Pauli matrices r J , \i % and a 1 . The r* act 
within the Dirac space of each two-component fermion, 
the a 1 act on SU(2) spin indices, and the fx 1 connect the 
two different nodes. The imaginary-time action can be 
written S = f cPxCe, with 



■ia u 



(8) 

where 7 M = (t 3 ,t 2 ,— t 1 ) for fi = 0,1,2, respectively, 



and * 



The observant reader will notice that 



we have dropped any explicit velocity anisotropy for the 
fermions; it is instead grouped with the other perturba- 
tions consistent with the microscopic symmetries repre- 
sented by the ellipsis. 

We defer consideration of the perturbations to the first 
two terms of Eq. JSJl until Sec. IIII Bl For now we simply 
drop them. It is immediately clear that the resulting the- 
ory has a much higher symmetry than that of the spin 
model. In addition to the t/(l)fl ux symmetry discussed 
above, there is an SU (4) flavor symmetry acting on the 
fermions. This symmetry is generated by the 4x4 trace- 
less, Hermitian matrices T a , where a — 1, . . . , 15. The 
T a can be expressed in terms of tensor products of the 
a 1 and // Pauli matrices; that is, they are linear com- 
binations of the basis {<r l , //, c 1 /^}. The action on the 
fermion fields is given by: 



* -> exp(iA a T a )$ 
§ -> vj/exp(-a a T a 



(9) 



Note that this is a flavor symmetry; that is, it rotates the 
4 two-component fermions into one another, but does not 
affect the Dirac structure. More precisely, [T a ,jfj] = 0. 
The remarkable consequences of this 5C/(4) symmetry 
are the main focus of this paper. 



B. The large-TV/ limit 

The field theory of Eq. JSJ) has a non-trivial confor- 
mally invariant fixed point that is not amenable to a di- 
rect analytical treatment. As with other such critical 
theories in 2 + 1 dimensions, the best that can be done is 
to deform the model to a limit where we do have control 



and use this to understand as much as possible about 
the case of physical interest. A useful and familiar anal- 
ogy is with the critical fixed point of the classical 0(n) 
model in three dimensions. To access this fixed point 
analytically it is necessary to study it in an e-expansion 
near four dimensions or in an expansion in 1/n directly 
in three dimensions. Both of these expansions similarly 
provide useful analytic access in the present problem as 
well. However, the fixed point that describes the theory 
in Eq. (JHJ) has no relevant perturbations [in contrast to 
the 0(n) critical fixed point]. Here we will we follow pre- 
vious wor ks§iLiiii2ii£i22i2ii and generalize the theory to a 
large number of fermion flavors by adding an extra in- 
dex to the Dirac field: *S? — > <5 Q , where a = 1, . . . , Nf/4. 
With this convention Nf is the number of two-component 
Dirac fermions and the flavor symmetry is enlarged to 
SU(Nf). Nf = 4 corresponds to the physical case of 
SU(2) spin. For Nf sufficiently large, it is reasonable to 
treat 1/Nf as a formal expansion parameter. Provided 
we take e 2 ~ 1/iV/, and in the absence of perturbations, 
the theory can be solved order-by-order in 1/Nf. This 
can be carried out simply in terms of diagrams, and is 
described in Appendix iBl 

It is believed that the large- Nf expansion describes a 
conformally invariant fixed point to all orders in 1/Nf 
This fixed point is the algebraic spin liquid. At Nf = oo 
the theory is scale-invariant and the fermions behave for 
most purposes as if they were free. (That the fermions 
are not truly free is apparent from the presence of op- 
erators that acquire an anomalous dimension even at 
Nf = oo; see the discussion of the gauge charge cur- 
rent in Sec. IIII Al This is analogous to the situation 
with the quadratic "mass" operator in the 0(n) model.) 
The 1/Nf corrections to this extreme limit correspond 
physically to incorporating gauge fluctuations, and one 
finds that various correlators acquire anomalous dimen- 
sions for Nf < oo. The usual justification for the pres- 
ence of a conformally invariant fixed point comes from a 
consideration of the above perturbation theory for an ar- 
bitrary correlation function. Due to the l/\q\ form of the 
photon propagator, the effective expansion parameter for 
this perturbation theory is easily seen to be dimension- 
less. Because there is no longer any scale in the problem 
(aside from a short-distance cutoff) it is natural to ex- 
pect that the large-Ay expansion describes a conformally 
invariant fixed point. Furthermore, a fermion mass can- 
not be generated perturbatively in 1 /Nf because all such 
terms break cither the SU(Nf) flavor symmetry, or par- 
ity and time reversal. 

Further insight is provided by a renormalization group 
(RG) approach perturbative in 1/Nf, which is also very 
useful as a tool for calculationiiiiSiSSiSiiSI One simply 
regards 1 /Nf as an exactly marginal perturbation to the 
Nf = oo fixed point, and calculates corrections to the 
properties of this fixed point as an asymptotic series in 
1/Nf. For technical purposes it is most convenient to 
implement a "field theory" RG, and first calculate some 
correlation function to the desired order in 1/Nf with a 
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fixed UV cutoff. Then we demand that this correlator 
satisfy the appropriate Callan-Symanzik equation, which 
is simply the mathematical statement that we can equiv- 
alently change the overall momentum scale k — > e~ e k, or 
rescale the fields. We can also include perturbations to 
the fixed point; then we must also rescale their coupling 
constants, and the resulting Callan-Symanzik equations 
allow us to calculate the flow equations for the couplings. 

In this language, the statement that the 1/Nf ex- 
pansion describes a scale-invariant fixed point to all or- 
ders can be put as follows: Set all perturbations to the 
fixed-point theory to zero. Then we can write a Callan- 
Symanzik equation for any correlator involving only the 
anomalous dimensions of fields. If these equations can 
all be satisfied order-by-order in 1/Nf, then the theory 
is indeed scale- invariant. While it has not been proven 
that this is the case, to our knowledge no inconsistency 
has been found. 



III. SYMMETRIES AND STABILITY OF THE 
SPIN LIQUID 

A. Symmetries and conserved currents 

The above considerations strongly suggest that for suf- 
ficiently large but finite Nf, somewhere in its parameter 
space the field theory has a conformally invariant fixed 
point smoothly connected to the Nf = oo fixed point. It 
is reasonable, although not certain, that this fixed point 
continues to exist for Nf — 4, and that a qualitative 
picture of its properties is provided by low-order calcu- 
lations in the 1/Nf expansion. This should be viewed as 
providing a definition of the algebraic spin liquid fixed 
point. 

The ASL then has the symmetry group 

•SGasl = (Conformal Symmetry) (10) 

X CxVxT X ST/^flavo, X C/(l)fl ux . 

This is a much larger symmetry than is present in the mi- 
croscopic model. In addition to the 5J7(4)fl avor , t/(l)fl ux , 
and conformal symmetries discussed above, we also have 
the discrete symmetries of charge conjugation (C), parity 
(■p), and time reversal (T). It is important to note that 
here we are referring to the symmetries of the continuum 
theory, and not, for example, to the time-reversal sym- 
metry of the spin model. The action of this operation 
on the field theory degrees of freedom will involve the 
continuum T combined with other operations. 

In Sec. IIII 51 we show that, in the large- Nf limit, all 
allowed perturbations are irrelevant and the staggered 
flux ASL is thus a stable phase. For Nf = 4 it is then 
likely that the sF state either remains stable or has only 
a small number of unstable directions. Since the physical 
case of Nf = 4 is of primary interest, we focus here on 
this case but use the large- Nf expansion to control our 
results. 



We wish to put special emphasis on the SU(4) flavor 
symmetry discussed in Sec. Ill Al which leads to a host 
of heretofore unnoticed consequences. In order to under- 
stand the relationship to the microscopic symmetries, it 
is useful to think in terms of the subgroup 

SU(2) spin x SU(2) nodal c SU(4), (11) 

where SU(2) sp i n is the physical spin and is generated by 
a 1 . S'L r (2) no dai consists of emergent symmetries and is 
generated by the Pauli matrices //. We refer to this as 
the "nodal" SU (2) because these flavor rotations involve 
the two distinct nodes arising from the staggered flux 
band structure and commute with SU(2) spin rotations. 
Certain discrete SU(2) noc i a x rotations are intimately tied 
to the microscopic lattice symmetries; this is apparent 
upon inspection of the symmetry transformation laws 
enumerated in Appendix lAl 

The conserved SU(4) flavor current is 

K = -i^T a «f. (12) 

This multiplet of operators transforms as an SU(4) ad- 
joint and a Lorentz vector. It is easy to show by explicit 
calculation in the Nf — oo theory that the two-point 
function of J?, falls off as 1 /x 4 , and J" thus has dimension 
2. Since conserved currents cannot acquire an anomalous 
dimension, this must hold for all Nf. 

The emergent Z7(l)fl ux symmetry also plays an impor- 
tant role. This symmetry is associated with the con- 
served gauge flux current 

j% = e^xd^ax- (13) 

This current also has scaling dimension 2 for all Nf. 
The conservation law dpjj* = is violated precisely 
by magnetic monopoles, so the emergent J7(l)fl ux sym- 
metry encapsulates the irrelevance of monopoles at low 
energies Monopole operators are those carrying a 
nonzero f7(l)fl ux charge. 

There is also a conserved U(l) gauge charge current 
G M = — i^j^ associated with the global gauge trans- 
formation "J — > . Since this gauge "symmetry" is not 
a true symmetry but rather the consequence of a redun- 
dancy in our choice of variables, it should not be surpris- 
ing that G^ is rather special. Making the infinitesimal 
change of variables a fl (x) — > a fl (x) + e M (x) in the func- 
tional integral leads via standard manipulations to the 
Maxwell equation, which we can regard as an operator 
identity: 

Gn = —ztu.v\du3\ + (more irrelevant terms). (14) 

As with all operator identities, the meaning of this equa- 
tion is that in any correlation function G^{x) can be re- 
placed by the right-hand side of Eq. i|14[l. as long as the 
other fields involved in the correlator are not close to x. 
This implies that G M has dimension 3, as can be verified 
by explicit calculation in the Nf = oo theory. Further- 
more, Gfj, should be thought of as a derivative of and 
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not a new truly independent operator. We note that this 
is a manifestation of the fact that the fermions are not 
quite free, even at Nf = oo. 



B. Stability of the spin liquid 

In order to assess the stability of the sF state at large 
Nf , we need only follow the argument of Ref . [l3 for the 
stability of the 7r-flux (ttF) state of an SU(N) magnet. 
The field theories in these two cases are identical, but 
different perturbations are allowed by the microscopic 
symmetries. As in Ref. ^2 we can group all operators at 
the sF algebraic spin liquid fixed point into two classes: 
those carrying t/(l)fl ux charge and those that do not. 
The first class is comprised of the monopole operators, 
which arc all strongly irrelevant at large Nf , with scaling 
dimensions proportional to Nf£* 

The monopole-free sector of the theory must be consid- 
ered in more detail, as it contains operators that would be 
relevant if allowed by symmetry. As in Ref. ll2l the poten- 
tially dangerous perturbations are fermion bilinears with 
zero and one derivative ( "mass" and "kinetic" terms, re- 
spectively). It is a simple exercise to show that all mass 
terms are forbidden by symmetry. In fact, this is true 
even if one only considers SU(2) spin rotations, time re- 
versal, and x- and y- translations. 

We must also consider the kinetic terms as these have 
dimension 3 + 0(1/Nf) and are exactly marginal at infi- 
nite Nf . Two such terms are allowed by symmetry. The 
first is simply the isotropic kinetic energy 



A" 



iai 



\d 2 



(15) 



This term has no effect, as it can be absorbed into the 
fixed point theory Eq. (JSJ) by rescaling the time coordi- 
nate. 

The second term is a velocity anisotropy for the 
fermions and cannot be removed by rescaling space and 
time. It is therefore important to know the l/Nf cor- 
rection to the scaling dimension of this operator, which 
takes the form 



K a 



tf#V + ioi) - r 2 (d 2 + iaa) (16) 



where it should be noted that the coefficient is 6. The 
RG flow of 5 can be calculated as a function of l/Nf and 
6. This was done by Vafek, Tesanovic, and Frana22iS^ 
to leading order in l/Nf and to all orders in <5; a more 
straightforward but equivalent calculation is discussed in 
Appendix [U] and reproduces the leading term of their 
result. We find 
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S + 0(S 2 /Nf,S/N 2 f ) 



(17) 



and the velocity anisotropy is therefore irrelevant for suf- 
ficiently large Nf. 



Other gauge-invariant operators in the monopole-free 
sector of the theory can be constructed by forming poly- 
nomials of the fermion fields and the gauge flux, and 
inserting covariant derivatives. It is easy to see that all 
such operators are trivially irrelevant in the Nf — > 00 
limit, and thus do not destabilize the algebraic spin liq- 
uid at sufficiently large Nf. 



IV. PHYSICAL OBSERVABLES 

In this section we examine two classes of field theory 
operators, focusing on their connection to observables in 
the spin model. The first class has slowly-varying corre- 
lations and gives rise to the competing orders within the 
algebraic spin liquid. The second class is comprised of 
the conserved currents J® and jff . 

Quite generally, any operator in the spin model is con- 
nected to the field theory by the relation 



O 



Spin Model 



7 , c i^Ficld Theory- 



(18) 



The meaning of this expression is that the long-distance 
correlations of the spin model operator are identical to 
those of the sum of field theory operators on the right- 
hand side. The Cj are nonuniversal coefficients, and 
generically Cj 7^ if and only if Op ic i d Theory transforms 
identically to the spin model operator under the micro- 
scopic symmetries - when this is the case we say the two 
operators are symmetry-equivalent. More precisely, all 
terms on both sides of Eq. I|18|) should transform in the 
same irreducible representation of the microscopic sym- 
metry group, and all field theory operators transforming 
in this representation will contribute. 



A. Fermion bilinears with enhanced correlations 

Here we shall be interested in the S£/(4) adjoint 

7V a = ~i^/T a ^, (19) 

as well as the scalar 



M 



(20) 



Rantner and Wen calculated the two-point function, and 
hence the scaling dimension, of one member of the N a 
multiplet to leading order in the l/Nf expansion^ - the 
operator they considered is symmetry-equivalent to the 
Neel vector. By ST/ (4) symmetry their result applies to 
the entire multiplet and its scaling dimension is 
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■0(1/N 2 f ) 



(21) 



The correlations of these operators are therefore en- 
hanced by gauge fluctuations - this is physically very 
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Field Theory 


Spin Model 


N\ , N\ 


( — l) ra: + 1 S r X S r + y , ( — l) Ty S r X S r + x 


N\ 




N b 


(-l) r »+ r »[(5i + 5s)(&-S4) 
+ (S 2 + S 4 )(Si-S 3 )l 


Nh , N 2 C 


( — l) ry S r ' Sr + y , ( — if" Sr ' ■S'r + a; 


N% 


[Si ■ (S 2 x S 4 ) - S 2 • (S 3 x Si) 
+ S 3 • (S 4 x S 2 ) -S 4 • (Si x S 3 )l 


M 


[Si ■ (S 2 x S 4 ) + S 2 ■ (5 S x Si) 
+ S 3 -(S 4 x S 2 ) + S 4 -(Si x S 3 )] 



TABLE I: List of observables in the spin model that are 
symmetry-equivalent to the N a and M fermion bilinears. For 
some of these we label the sites around the plaquette with 
lower- left corner at r by the numbers 1, . . . , 4. Precisely, Si = 

S r , S 2 = S r +x ; S 3 = S r + x +y and S 4 = S r +y • 



reasonable, since the gauge force tends to bind the oppo- 
sitely charged 'J and particles. 

Although they are not related by any obvious symme- 
try, M and N a have the same scaling dimension to all 
orders in 1/Nf, this is shown in Appendix IDl At present 
it is not clear if there is a deeper reason (e.g. some hid- 
den symmetry) for this equality. Therefore, while this 
statement may hold at finite Nf, it may also be merely 
an accident of the large- Nf expansion that is destroyed 
by effects nonperturbative in 1/Nf. Whether or not M 
and N a have the same scaling dimension, the large- Nf 
expansion indicates at least that the correlations of M 
are also enhanced by gauge fluctuations. 

At large Nf, M and A^ a are the most relevant operators 
and give the dominant long-distance correlations. They 
are also in some sense the most natural instabilities of 
the algebraic spin liquid, although a proper treatment of 
this issue requires consideration of the parameter space 
around the hxed point and not only the ASL itself. 

It is useful to find physical observables symmetry- 
equivalent to N a and M as their correlations may decay 
slowly enough to be readily observable. This exercise is 
easily carried out making use of the transformation laws 
of Appendix ^ here we summarize and discuss the re- 
sults. It is convenient to group the N a into three classes 
depending on their matrix structure: 

7V^ = -if/iV*, (22) 
Nb = -i*cr*, (23) 
N l c = (24) 

Symmetry-equivalent spin operators are listed in TableQ] 
These operators are order parameters for several slowly 
fluctuating competing orders. It is remarkable that these 
apparently unrelated fluctuations are perfectly balanced 
within the algebraic spin liquid. Two of these observ- 
ables are quite familiar: N^ is the Neel vector, and 
$vbs = {Nc, N*) is the valence-bond solid order param- 
eter. When $vbs has an expectation value along the x 



or y axis the columnar VBS state (Fig.^3) results, while 
if it points at 45° from the axes the resulting state is the 
box VBS (Fig.njU). We note that this unification of Neel 
and VBS order is quite different from the situation at the 
recently elucidated quantum critical point between Neel 
and VBS states^SiS where the order parameters have 
different scaling dimensions [in the case with full SU (2) 
spin symmetry] and are not related by symmetry. 

The other order parameters are rather unusual and cor- 
respond to more exotic ordered states. Nb transforms 
like a plaquette-centered spin at q = (tt, it). The opera- 
tors iVV form the order parameter for a kind of triplet 
valence bond solid that breaks spin rotations but not time 
reversal. 

Nq is somewhat more familiar. If we denote the slowly- 
varying Neel field by the unit vector field n, we can define 
the density of the familiar topological Skyrmion configu- 
rations by writing 

1 ,dn <9ro N 
Pfl = 4^ n ' ( fl£ X fl*>- (25) 

As usual, if we take periodic boundary conditions in 
space, J cPrps is an integer that changes from one time 
slice to the next upon encountering an instanton defect 
in the Neel field (a hedgehog). Now, A^^, and the corre- 
sponding symmetry-equivalent spin operator transform 
like the q — (tt,tt) component of ps- 

Finally we turn to M, which is odd under time rever- 
sal and reflections, and does not transform under any of 
the other microscopic symmetries; it is thus symmetry- 
equivalent to the uniform component of the scalar spin 
chirality. Furthermore, under the symmetries of the ASL 
fixed point, M transforms exactly like a Chern-Simons 
term for the gauge field a M , and if M is added to the La- 
grangian (via spontaneous symmetry breaking), a Chern- 
Simons term will also be generated. The resulting state 
is a chiral spin liquid supporting gapped spinons with 
fractional statisticsi 32 i 33 i 34 



B. Observables for the conserved currents 

Now we shall discuss symmetry-equivalent spin oper- 
ators for the two conserved currents J° and iff. This 

fj, j fi 

is most interesting for the information obtained about 
the general structure of the sF state. As discussed in 
Sec. IIII A1 the corresponding correlations decay as 1/r 4 
and may be rather difficult to observe. Furthermore, we 
will see below that a given microscopic operator may be 
symmetry-equivalent both to a member of one of the cur- 
rent multiplets and another, more relevant operator. In 
order to predict that a particular correlation function of 
the microscopic model decays as 1/r 4 we need to know 
that this does not happen. For numerical simulations and 
experiments, it is undoubtedly better to begin by looking 
for the stronger correlations discussed above, and to con- 
sider the "fine structure" of the conserved currents only 
as a second step. 
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Rather than systematically considering every compo- 
nent of the conserved currents, we only highlight some of 
the most interesting cases. We begin with the gauge flux 
current . In the (continuum) Hamiltonian its compo- 
nents correspond to the magnetic flux <1>b = Jq and the 
electric field = ieijjf ■ We find that the electric field is 
symmetry-equivalent to a staggered scalar spin chirality, 
taken along lines of three adjacent lattice sites: 

(-l)^+r v ) Sr _ x . ( Sr x Sr+x) „ (26) 

-(-l)( r '+ r ')S r _ 9 .(S,xS r+9 ) ~ E y . (27) 

Labeling as above the four sites of the plaquette with 
lower-left corner r by the numbers 1, . . . , 4, two operators 
symmetry-equivalent to the magnetic flux are 

(-l)( p -+ r ») ( (S 1 -S 2 )(S 3 -S 4 ) 

+ (S 2 -S 3 )(S 1 -S 4 )) (28) 

and 

(-ljfr-^jSj • S 3 )(S 2 ■ S 4 ) ~ * B . (29) 

Subtracting Eqs. I|28[l and l|29|) we also see that the 
q = (tt, tt) component of the usual ring exchange op- 
erator (i.e. that obtained from the Hubbard model at 
order t 4 /U 3 ) is symmetry-equivalent to $b- 

Now we move on to the 5C/(4) flavor current J". The 
spin at q — (it, 0) and q = (0, tt) is symmetry-equivalent 
to two spatial components of the current: 

{-l) r *S r ~ *V + t 2 )/zV*, (30) 
{-l) r *S r ~ * t (-r 1 +r 2 ) M 1 <T*. (31) 

Rantner and Wen calculated the leading 1/Nt correc- 
tions to the correlations of these quantities and found 
no anomalous dimension. 11 This result is explained by 
SU(A) symmetry (and conformal invariance), which im- 
plies these operators have dimension 2 to all orders in 
1/Nf (see Sec. lmAl 

The two components of the VBS order, which already 
made an appearance above in the JV a multiplet, are also 
symmetry-equivalent to two of the SU (4) conserved den- 
sities: 

(-iy*s r - s r+x ~ *y*, (32) 

(-1)^ S r ■ Sr+y ~ *V*- (33) 

Note that there is no inconsistency in the fact that the 
VBS order appears in two distinct multiplets of the field 
theory. This simply means that both field theory opera- 
tors contribute to its long-distance correlations; that is, 

{-l) r *S r ■ S r+X ~ d-i-VV * + c 2 *V * + ■ ■ • , (34) 

where ci : 2 are nonuniversal constants as in Eq. I|18f) . This 
is an example where it is clearly not true that VBS cor- 
relations fall off as 1/r 4 simply because the VBS order 



parameter appears in J" This possibility must be con- 
templated for other operators symmetry-equivalent to 
conserved currents, and in general it is necessary to con- 
sider symmetry-equivalent field theory operators beyond 
the fermion bilincars. In particular, monopole operators 
carry nontrivial quantum numbers and may play into 
these considerations 

The Skyrmion density ps is symmetry-equivalent to 
one of the conserved densities: 

p s ~*V*- (35) 

It is very interesting that ps is also conserved at the de- 
confined critical point between the Neel and VBS states, 
where it corresponds to the magnetic flux of an emergent 
gauge field. We let be the gauge flux current as de- 
fined in Ref |23- There it is denoted as jj*; here we call 
it Sfj, to emphasize that it is quite distinct from the j~* 
defined in Eq. (|13f) - in particular, the two currents are 
not symmetry-equivalent. We do find, however, that 5^ 
is symmetry-equivalent to the following components of 
the SU(4) flavor current: 

50 - *V 3 *, (36) 

51 ~ -J=*t (r l_ r 2 )/i 3^ ) (37) 

5 2 ~ +r 2 )pH. (38) 

v2 

The presence of a 45° rotation in the r-matrix structure 
is to be expected, since the continuum coordinates of the 
sF state are rotated from the lattice axes by 45° (see 
Appendix This is not the case for the continuum 
theory of Ref.^ 

Remarkably, then, the conserved Skyrmion current 
that plays such a key role at the Neel- VBS critical point 
is also a conserved current of the sF algebraic spin liquid. 
Furthermore, it is contained within the larger structure 
of the SU{A) flavor symmetry It would be interesting 
to see if there is a natural route between these two fixed 
points - the U(1)r uk symmetry would have to be bro- 
ken by monopole proliferation, and the SU (4)fl avor would 
need to be broken down to SU(2) sp i n x U(l), where the 
U(l) corresponds to Skyrmion number conservation at 
low energies. On an even more speculative note, per- 
haps other interesting fixed points, so far undiscovered, 
also have the seeds of their structure hidden within the 
algebraic spin liquid. 



V. CONSEQUENCES FOR THE tt-FLUX STATE 
OF AN SU(4) HEISENBERG MODEL 

Recently, Assaad has carried out a quantum Monte 
Carlo study of an SU (4) Heisenberg antiferromagnet on 
the square latticed The results available to date are con- 
sistent with observation of the 7r-flux algebraic spin liq- 
uid, first studied in the large-A^ limit by Affleck and 
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Marstoni&i This state has a very similar structure to the 
sF ASL considered up to now in this paper - the pri- 
mary difference is that now we have a microscopic SU (4) 
spin symmetry, and the microscopic lattice symmetries 
act differently on the continuum Dirac fields. Here there 
is an emergent SU (8) flavor symmetry, and there are 64 
fermion bilincars which have correlations enhanced by 
gauge fluctuations. One of these is the SU(4) analog of 
the Neel vector - Ref. found that its correlations fall 
off very slowly, as l/r a , where a « 1.1 — 1.2. 

The SU (8) symmetry allows us to make the highly non- 
trivial prediction that certain other observables should 
have the same long-distance correlations. This will hold 
in the irF state provided that the most relevant operators 
for Nf = 8 are indeed the fermion bilinears discussed be- 
low, as is suggested by the large- Nf expansion. If some 
other multiplet of operators is more relevant, a similar 
set of predictions will hold, but for different observables. 
It is particularly important to consider monopole opera- 
tors in this context, since they carry nontrivial quantum 
numbers** and may have relatively low scaling dimen- 
sion; this issue will be considered in more detail in a 
forthcoming paper;*- With this one caveat in mind, the 
results of this paper can be tested numerically, and it 
should thus be possible to determine rather conclusively 
whether or not the irF state has indeed been observed in 
the model of Ref. 0. 

The model is defined in terms of the slave fermions f ra , 
where a — 1, ... ,4, and we choose the local constraint 
fra fra = 2- Here the SU(4) spin rotations are a micro- 
scopic symmetry, again generated by the 4x4 matrices 
T a , with a = 1, . . . , 15. The action of an SU(4) spin ro- 
tation on the fermions is f ra — > [exp(i\ a T a )] a ^f rl 3. We 
define the Hermitian SU (4) spin operator 

&r = fraTapfrP' 0*9) 

and the Hamiltonian takes the form 

%(4)=^S r X, (40) 

{rr' ) 

with J > 0. 

The quartic spin-spin interaction can be decoupled 
with a compact U(l) gauge field, and one can consider a 
mean- field theory (exact at infinite Nf) where this gauge 
field becomes a non-fluctuating classical background. In 
the mean-field 7rF state there is a gauge flux of ir through 
every plaquette, and the mean-field Hamiltonian has the 
same form as Eq. (jSJ) with t = A. The discussion now 
proceeds almost identically to that for the staggered-flux 
state above. The universal physics of the 7rF state, in- 
cluding fluctuations, is encapsulated in the lattice gauge 
theory Hamiltonian 

Hkf = h E elr> ~ K ^2 cos ( curla ) ( 41 ) 

{rr') D 

-*E E W + i-^'^fLe^-'fr^ + a-c.}. 

r£A r'n.n.r 



The low-energy effective field theory can be described by 
eight massless two-component Dirac fermions minimally 
coupled to a noncompact £7(1) gauge field. The field 
theory is thus the same as for the sF state, except that 
now Nf = 8. 

There are some differences between the ttF and sF 
states. First of all, the microscopic symmetries of the 
lattice model act rather differently on the continuum 
Dirac fields - for the 7rF state, the symmetries are enu- 
merated in Appendix A of Ref. 12. Also, the SU(2) 
spin version of the 7rF state is not a U(l) spin liquid 
at all, but instead has a gapless SU{2) gauge bosonii* 
Finally, the model Eq. I|40(l has an additional discrete 
global "charge-conjugation" symmetry (called C) with no 
analog in SU(2) spin models. C is defined as a particle- 
hole transformation of the spinon operators: f ra — > f} a . 
[Note that in an SU (2) spin model this is equivalent to a 
particular SU(2) spin rotation and is thus not a distinct 
symmetry.] 

The 7rF state has an SU (8) flavor symmetry, and there 
is again a useful decomposition into the subgroup 

St/(4) spin x SU(2) nodal c S77(8) flavor . (42) 

The spin SU(4) is generated by the T°, and the nodal 
SU(2) is generated by the fi z Pauli matrices. We define 
the generators of SU(8) to be T A , where A = 1, . . . , 63. 
The T A can be expressed as tensor products of the T a 
and /Li*. Proceeding as in the sF state, we define the 16- 
component fermion field 5*, with all the flavor and Dirac 
indices sup pressed. [The notation here is identical to that 
of Ref. Il2l except for the very minor difference that here 
the SU (4) spin index is suppressed and \& is written in- 
stead of ^> a -] The action of the SU(8) symmetry is then 
— > exp(i\ A T A )ty . As above, we can form the 64 bilin- 
ears with correlations enhanced by gauge fluctuations: 

N A = -iVT A q, (43) 
M = -if*. (44) 

As for the sF state, it it convenient to break the N A 
operators into three classes: 

N a / = -i^T a f/^f, (45) 
Ng = -i*T Q *, (46) 
N l c = -if/i 1 *. (47) 

The goal here is to find microscopic operators with the 
same transformation properties as the continuum bilin- 
ears under the microscopic symmetries - this is easily ac- 
complished by making use of results in Ref. It should 
then be possible to numerically measure the correlations 
of these observables. 

At this point, it would be natural to find spin operators 
built from S% that transform as the various components 
of N A and M. However, this is not the most convenient 
way to proceed, since the simulation of Ref. works 
directly in terms of the slave fermions and the micro- 
scopic gauge field used to decouple their quartic inter- 
action. If we restrict our attention to spin operators, 
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Field Theory 


Lattice 




(-!)"" [eX P ((-iy* +r ^)W r , r + yfLT^f r+y 

+ H. c] 


N a / 


(-l) r -[exp )W r ,r +x fLT^f r+Xi 
+ H. c.j 


N a / 


(-i) (rx+ry) fLT^f rd 


N% 


(-l) {rx+ry) [(W h3 fLTS f 3 + H. c.) 

+ {W 2A fl a T^f ig + H. c.)] 


Nh 


(-iyy [exp ((-l)^3jr )W r , r+ yfLfr+v, a 

+ H. c.] 


N 2 C 


c-ir [exp ((-ir + ^f )Wr,^i/ r+Ii0 

, TT 1 

+ H. c. 


Nh 


('_-n(>-x+>- H ) /-t /• 

V A 7 J raj r a 


M 


(-l)< r - +P »J[(Wl, S /L/sa+ H.C.) 

+ (W / 2,4/ 2 t Q / 4Q + H. C.)] 



we denote a rr / . The two gauge fields are related by 

a r,(r+a;) = 5r,(r+a>) + (-l) rx+rH ^ — ^ , 

«r,(r+y) = flr,(r+y)+ (-1)^+''" (^) . (48) 

In terms of a rr ', the hopping term in the Hamiltonian 
Eq. 141(1 becomes 

-V2tJ2 (fLe- l ~ a -' f r , a + H. c.) . (49) 

{rr') 

The factor W r y in Table [H] is a function of the vec- 
tor potential included to keep the lattice bilinears gauge- 
invariant. For r and r' nearest neighbors we simply have 
the usual exponential 

W r y = cxp(— ia r . r >). (50) 



TABLE II: Continuum fermion bilinears with enhanced cor- 
relations and their lattice counterparts for the 7r-nux state 
of the SU(4) Heisenberg model discussed in the text. The 
factors W r r i encode the dependence on the lattice vector po- 
tential and are defined in the text. As in Table [I] for Ng and 
M we consider the plaquette with lower-left corner at r and 
label the sites around it with the numbers 1, . . . , 4, proceeding 
counterclockwise from r. 



most of the resulting observables are products of two or 
more S" and are therefore of quartic or higher order in 
the fermions. This is undesirable, most simply because 
such operators are rather difficult to deal with numer- 
ically. Furthermore, it is conceivable (because Nj = 8 
may be rather large) that the microscopic slave fermions 
are good variables and rather accurately represent the 
long-wavelength degrees of freedom. If this is the case, an 
operator quartic in f ra will have very little overlap with 
the continuum bilinears N A or M - precisely, the coeffi- 
cient of these operators in Eq. I|18|) will be dominated by 
that of an appropriate four-fermion term. So while the 
bilinears should indeed give the dominant long-distance 
correlations, it may be necessary to go to unreasonably 
large distances to overcome the small prefactor. 

To avoid these problems, we instead consider bilinears 
of the lattice slave fermions 4£ In most cases these involve 
products of fermions on different lattice sites, so in or- 
der to write gauge invariant operators it is necessary to 
include an appropriate dependence on the vector poten- 
tial. In Table HT1 we enumerate lattice bilinears and their 
continuum counterparts - in addition to transforming 
identically under the microscopic symmetries, the lattice 
observables also have the desirable property that they 
reduce exactly to the corresponding continuum operator 
upon taking the naive continuum limit. It is important to 
note that the form of these operators depends on the pres- 
ence of the explicit background flux in the Hamiltonian, 
Eq. (|41|l . This background is not present in Ref . [ijj this 
is simply due to a different choice of gauge field, which 



We also need to define W for r and r' next-nearest neigh- 
bors. In that case, to retain as much symmetry as possi- 
ble we must sum over the two shortest paths connecting 
r and r'. That is, if r' = r + x + y, we define 

W r y = exp ( - ia r ^ r+x ) - ia^ r+x - ) ^ r+x+y ~ ) ) (51) 

+ ex P ( ~ i- a r,(r+y) — ^ a (r+y) ,(r+x+y)) ■ 

While it is not directly useful for obtaining numerical 
results, it is interesting to discuss the orders in the 7rF 
state in more physical terms. It is easily seen that iV^' 3 is 
the SU(A) version of the Neel vector. Also, N^ 2 together 
form the order parameter for the columnar and box VBS 
states, as in the sF case. Specifically one finds 

{~iyyS a r S a r+y ~ Nh, (52) 
{-iy*S a r S a r+x ~ N 2 C . (53) 

Nq is an order parameter for a kind of C-breaking 
state that has been studied for a large class of SU(N) 
magnets. 36 It should be noted that Nq breaks time re- 
versal as defined in Ref. (there time reversal sends 
S a — > — S a ), but does not in the conventions of Ref. l3rl 
Finally, as in the sF case, M transforms identically to 
a Chern-Simons term for the gauge field, and if sponta- 
neously generated will lead to a chiral spin liquid. 



VI. CONSEQUENCES FOR UNDERDOPED 
CUPRATES 

In our view, one of the more promising routes toward 
a theoretical understanding of the underdoped cuprates 
views the pseudogap regime as a doped spin-liquid Mott 
insulator. This line of thinking began with the ideas of 
Anderson— and has subsequently been developed by many 
othersi 37 i 38 

The current state of these ideas has been discussed 
recently by Senthil and Lee£ For our present purposes, 
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FIG. 2: Schematic temperature-chemical potential phase di- 
agram of an underdoped cuprate superconductor. For /i < fi c 
the holes are gapped and the ground state is a spin liquid 
Mott insulator, while when fi > /i c a finite density of holes 
has entered the system and the ground state is a d-wave 
superconductor. The regime labeled "QC" is controlled by 
the Mott quantum critical point separating these two phases, 
while the FS regime is best thought of in terms of fluctu- 
ating d-wave superconductivity. The dashed line represents 
the finite-temperature behavior of a real underdoped material 
(with a superconducting ground state) at fixed doping. 



the key point of this picture is the close proximity of 
the underdoped d-wave superconductor to a Mott tran- 
sition to a spin liquid insulator. This is most clearly 
understood by thinking about the phase diagram as a 
function of hole chemical potential and temperature, as 
shown in Fig. (from Ref. Consider an underdoped 
material in the superconducting state, and imagine rais- 
ing the temperature - this is represented by the dashed 
line in Fig. Above T c there is a region that can be 
described as a phase-fluctuating d-wave superconductor 
(labeled by "FS" in Fig. |2J, and at still higher temper- 
atures the physics is controlled by the quantum critical 
point between the spin liquid and the superconductor. 
We shall be interested in this "high-temperature pseudo- 
gap" region. 

We take the Mott insulator in question to be the stag- 
gered flux algebraic spin liquid^Si^i Furthermore, we 
consider the simplest scenario where this state is a stable 
phase. In this case, the Mott insulating part of the phase 
diagram is controlled by a critical state with a dynamic 
critical exponent z = 1. The Mott transition is described 
by the condensation of a doubly-degenerate parabolic 
band of charge e bosonic holons. An important point is 
that the charge sector essentially has z — 2, and at long 
distances the holon velocity goes to zero while the spinon 
and photon velocities go to a constant. This means that 
the holons move very slowly and couple only weakly to 
the spin sector without strongly influencing it. Therefore 
the spin sector of the quantum critical regime should not 
change substantially from the finite-temperature physics 
of the algebraic spin liquid. It is very important to note 
that this conclusion will hold only at temperatures above 
the crossover to the fluctuating superconductor (or "FS" ) 
regime. Below this temperature, the coupling between 



the spin and charge sectors will be very important. 

It should then be possible to test whether the staggered 
flux algebraic spin liquid is relevant to the cuprates by 
probing the high-temperature pseudogap for signatures 
of the ASL fixed point. What are the upper and lower 
temperature scales defining this region? Physically, one 
needs a large window of length scales where the spin sec- 
tor is coherent and controlled by the ASL fixed point, but 
the charge sector is rather incoherent and sufficiently far 
from condensing that it does not substantially influence 
the spin sector. Presumably the upper temperature is 
the pseudogap scale T* , below which the Knight shift de- 
creases, signaling the formation of spin singlets. Because 
the uniform magnetization is a conserved density (one of 
the J"), it can be seen by a standard scaling argument 
that the Knight shift in the ASL is expected to decrease 
linearly with decreasing temperature^ The lower tem- 
perature scale is presumably bounded below by T v , the 
onset temperature for the Nernst signal 4 2a£ The onset 
of the Nernst signal can probably be identified with the 
onset of a substantial degree of phase coherence in the 
charge sector. 

To take a specific example, in YBa 2 Cu306+2; T v is 
about 110 K for x = 0.5^i while, for an x — 0.53 sample, 
the Knight shift increases linearly from this temperature 
up to 300 K (the highest temperature measured)^ So 
in this material it is reasonable to look for algebraic spin 
liquid physics at least in the range 110K < T < 300K, 
and perhaps at even higher temperature. This window is 
expected to grow with underdoping, which increases T* 
and decreases T v ^ 

Within this temperature range the properties of the 
algebraic spin liquid can be probed by looking at the 
dynamic correlations for frequencies uj <C cJ. The en- 
ergy cJ, where J is the exchange energy and c is a num- 
ber of order unity, plays the role of a high-energy cutoff 
above which the physics is presumably nonuniversal. If 
one further restricts to frequencies uj 3> T, it is possi- 
ble to probe the zero-temperature critical ground state 
and avoid complicated issues of critical dynamics. Fur- 
thermore, various quantities will exhibit critical scaling 
for all frequencies uj <C cJ. It is likely that the sim- 
plest test of this physics would be to look for scaling in 
the Q = {it, 7t) magnetic neutron scattering. Specifically 
the imaginary part of the dynamic spin susceptibility at 
(■7T, it) is expected to satisfy 



x"(q,",T) 



1 



{{q-QY+uj*) 2 



2-7) f 



T 



Q\ uj_ 

' T 



The exponent r\ is not known. A rough guess may be 
obtained from studies of variational wavefunctions for the 
sF spin liquid. In particular the well-studied projected 
nearest-neighbor d-wave BCS state might be expected 
to capture the physics of the sF spin liquid. From the 
known resultpi^^ on the equal-time spin correlations in 
that wavefunction one extracts rj sa 0.5. The scaling form 
above also has direct implications for the NMR relaxation 
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rate 1/Ti at the Cu site [which is sensitive to the (%, ii) 
spin correlations]. We have 1/Ti ~ T* 7 with 77 roughly 
about 0.5. 

If this type of scaling is seen, it would be important 
to think about whether the SU (4) symmetry can be ex- 
plicitly tested by probing the other observables with en- 
hanced correlations. In particular the power law VBS 
correlations can possibly be looked for. 

Very recently, scaling in the Q — (tt, tt) magnetic neu- 
tron scattering has been observed in underdoped sam- 
ples of YBa 2 Cu 3 06+x, for both x = 0.5 (T c = 59K) and 
x = 0.35 (T c = 18K)M In both samples cu/T scaling 
is seen in x"(q, uj,T)/x"(q,u 7 T = 0). Further analysis 
of the data is needed to understand whether this scaling 
may be related to algebraic spin liquid physics. 

We note that scaling has also been reported in very 
lightly doped La2-a;Sr x Cu04 fRefs. l44l andl45i) and more 
recently in Li-doped La2Cu04 (Ref . l46|) . Specifically the 
q-integrated scattering intensity (which is dominated by 
the signal near Q) shows u>/T scaling. However, the pref- 
actor decreases with increasing frequency unlike that ex- 
pected from the scaling form abovei^ Further, there is 
no sign of scaling in the momentum dependence of the 
scattering. We therefore think it unlikely that the ui/T 
scaling reported in Refs. l4H45Tl4^ is due to any underlying 
algebraic spin liquid. However, it would be very interest- 
ing to look for scaling in moderately doped samples in the 
high-temperature pseudogap regime (as op pos ed to the 
very lightly doped samples studied in Refs.0 and l45T) . 

It is important to note that, even in this relatively 
simple picture where the algebraic spin liquid is stable, 
the critical scaling may be modified by the presence of 
weakly irrelevant perturbations. Particularly worrisome 
is the fermion velocity anisotropy, which is known to be 
rather large in the superconductor. However, it is not 
at all clear that the low-temperature anisotropy in the 
superconductor is to be identified with the anisotropy in 
the spin sector at high temperatures. It may well be 
the case that the spin sector flows close to the isotropic 
ASL fixed point, but at lower temperatures the charge 
sector could induce the anisotropy that obtains in the 
ground state. Similar issues arise regarding the location 
of the nodal points in momentum space - in the sF state 
these arc fixed at q = (tt/2,tt/2), but this is not the 
case in the superconducting state. While these issues 
remain somewhat mysterious at present, we would like 
to emphasize that the whole approach of thinking about 
the underdoped cuprates in terms of an algebraic spin 
liquid is only useful to the extent that one comes near 
the fixed point, which is isotropic. 



VII. DISCUSSION 

In the analysis of this paper, we have assumed that 
the most relevant operators in the sF and irF spin liq- 
uids are the iV° and M fermion bilinears; these opera- 
tors therefore give rise to the dominant power-law cor- 



relations. While this scenario will certainly hold for suf- 
ficiently large Nf, it is not known whether it continues 
to hold for the interesting cases of Nf — 4 or 8. In par- 
ticular, monopole operators carry nontrivial flavor (and 
other) quantum numberS(2i^ and if some such multi- 
plet of operators becomes more relevant than, say, iV a , 
it will dominate the long-distance correlations. We note 
that if this does in fact happen, the main features of the 
results discussed here will still hold. Specifically, there 
will be a set of superficially unrelated and slowly- varying 
competing orders all exhibiting the same power-law de- 
cay and unified the by the emergent SU(Nf) symmetry. 
The observables involved will, however, be different from 
those discussed here. Furthermore, even if N a and M do 
not give rise to the dominant competing orders, they will 
still give rise to slowly-decaying correlations. It should 
be noted that these issues do not affect the prediction of 
lu/T scaling in the q = magnetic scattering. 

Many aspects of the staggered flux spin liquid state are 
very reminiscent of the physics of the more familiar one- 
dimensional critical spin liquids. The most striking simi- 
larity is perhaps in the criticality itself; indeed, the stag- 
gered flux and other algebraic spin liquids can be stable 
critical phases in two dimensions, much like their better- 
known one-dimensional counterparts. In both cases the 
spin correlations are described by nontrivial power laws 
with large anomalous dimensions. More technically, as 
discussed in detail in this paper, the sF spin liquid is 
conveniently analyzed in terms of fermionic S = 1/2 
spinon variables. Similarly, a fermionic description is of- 
ten a useful technical device in analyzing the physics of 
one-dimensional spin liquids. Finally, there is similar- 
ity in how the semiclassical instantons are represented 
in terms of the fermions. In d = 1, for instance in the 
antiferromagnetic XXZ S = 1/2 model, the semiclassi- 
cal instantons are just 2ir phase slips. In a fermionic 
representation (obtained via Jordan- Wigner transforma- 
tion) these are umklapp processes where a right mover 
becomes a left mover or vice versa. In the d — 2 spin 
model, the semiclassical instantons change the Skyrmion 
number associated with the Neel vector configuration. In 
the staggered flux spin liquid these are again described as 
operators that move a fermionic spinon from one node to 
the other (the northeast movers to the northwest movers, 
for instance). Perhaps these similarities can be exploited 
toward deepening our understanding of such nontrivial 
two-dimensional algebraic spin liquids. 

We note that the spin-charge-separated variables used 
to describe the algebraic spin liquids considered here have 
no obvious a priori connection to the pattern of compet- 
ing orders arising within these states, and it is remarkable 
that they lead to this kind of physics. In fact, the ob- 
servables with slowly-varying fluctuations correspond to 
bilinears of the fermions. It is not known how to formu- 
late a field theory for the sF state where these variables 
are in some sense the fundamental fields. This points 
out that, in doing phenomenological modeling of strongly 
correlated systems, one should be cautious about simply 
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introducing new fields by hand for the slowly fluctuat- 
ing observables - these are not necessarily the variables 
that will naturally lead to a correct description of the 
underlying physics. 

Many theoretical issues remain to be addressed if a 
solid connection is to be made between the picture ad- 
vocated here and experiments in the cuprates. We feel 
the most serious of these involve coupling to the charge 
sector, which so far has not been carefully taken into ac- 
count. Specifically, it will be important to understand 
the physics of the doped algebraic spin liquid at all tem- 
peratures, not only in the range where the spin sector 
should be controlled by the undoped fixed point. The 
zero-temperature fate of a doped ASL is a question that 
also merits further exploration. The most common view 
has been that <i-wave superconductivity obtains imme- 
diately at T — upon introduction of a finite density 
of charge carriers. It is particularly intriguing to ask 
whether some of the exotic character of the ASL can sur- 
vive down to T — even in the presence of doped holes, 
possibly leading to exotic metallic states. 
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APPENDIX A: CONTINUUM FIELDS AND 
MICROSCOPIC SYMMETRIES 



Here we provide a discussion of the continuum limit 
of the staggered-flux mean-field state, and the resulting 
action of the microscopic symmetries on the continuum 
fields. The same procedure was discussed in Appendix 
A of Ref. [H for the SU(N) tt-Aux state. Much of the 
analysis is identical, but for completeness we reproduce 
it here. It is important to note that the final results 
differ because the symmetries in the 7rF and sF states 
act differently on the lattice spinons and hence also on 
the continuum fields. 

The starting point is the mean-field Hamiltonian of 
Eq. JSJ. We choose the four-site unit cell labeled by 



where 



v, = 



0. 


i 


= 1, 


X, 


i 


-2, 


x -i 


y, i 


= 3, 


y- 


i 


= 4. 



(Al) 



The spinon operator at site (R, i) is denoted / ma- 
It is a trivial exercise to go to momentum space and 
solve Eq. JSJ; in the reduced Brillouin zone k x ,k y £ [0, ir) 
one finds gapless Fermi points at Qo = (ir/2,ir/2). Near 
this point the dispersion can be described by 4 two- 
component Dirac fermions. It is convenient to denote 
these by tp^ a (R). Here a = 1,2 and a = 1,2 are the 
ST/(4) flavor indices [a is simply the 5/7(2) spin index]. 
Also, A = 1, 2 labels the two components of each spinor 
(this is usually suppressed). These fields are related to 
the lattice spinons as follows: 







i>l a (R) - 




i>L{R) ~ 


_ e -w 


rtJR) - 


_ e -w 

2^21 



lQo R (f R 2a - faux), (A3) 
-^(/ata + W, (A4) 



-e lQo R (fma~ fR3a), (A5) 



where I is the lattice spacing. 

We can now set t = A to remove the velocity 
anisotropy, and put it back in as a perturbation as dis- 
cussed in Sec. IIII Bl Note that the above results do not 
depend on t/ A and are identical to the case of the 7r-flux 
state. In momentum space the continuum Hamiltonian 
takes the form 



d 2 q 
(2tT 2 ) 



V>L(q)(?i^ + i2T 2 )4> aa (q), (A6) 



where we have chosen units to set the velocity to unity, 
and t % are the usual Pauli matrices acting in the two- 
component Dirac "spin" space. Here we use the following 
rotated coordinates: 

5i = -^(ix + qy), (A7) 
1 / 

92 = -y={-q x + q y ). 

It is convenient to work with the eight-component ob- 
ject 



* = 



^12 
^21 
V "022 / 



(A8) 



(R,i), with R = 2n x x + 2n y y and r(R,i) 



R 



The generators of flavor SU(4) can be expressed as ten- 
sor products of SU(2) sp i D and SU(2) DO d a i generators (see 
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Sec. IIII A"|) . Using the convention specified in Eq 
this can be expressed in matrix notation; for example, 



(A9) 



We now quote the action of the microscopic symme- 
tries on the lattice and continuum fields, including the 
spinons as well as the magnetic flux and electric field. 
To simplify the form of the results, we often make an 



additional global gauge transformation f r 



i<t> 



fra in 



going from the lattice to the continuum transformation 
laws. The lattice symmetries discussed below generate 
the full space group of the square lattice. 

x -translations. Translations by one site in the x direc- 
tion act on the lattice spinons as follows: 



fra > £r(io~ )apfl +x jj, 
fra ~* tr(io- 2 ) a pf r + x ,l3, 



(A10) 



where 



I +1 r e A, 
1 -1 r S B. 

The resulting continuum transformation law is 

* -» [*t (ir 1 )(ia 3 )] T ', 
* f -» [(ia 2 ){iT 1 )V] T . 



(All) 



(A12) 



The electric field and magnetic flux both change sign un- 
der translation by one lattice site. 

Rotations. We choose to make a ir/2 counterclock- 
wise rotation about the point (x + y)/2, which lies 
at a plaquette center. Under this operation we have 
r — > r' = (— r y + 1, r x ), and the action on the spinons is 



fra * ^rfr'a- 



(A13) 



In the continuum 



*(*) -> exp (| (^)) exp (^)*(R')- (A14) 

Under this operation, the electric field is a vector, and 
the magnetic flux is a scalar. 

Reflections. We consider the reflection r — > r' = 
(— T x ,Ty). The spinons transform trivially: 



/ra * fr'ai 

resulting in the continuum expression 



(A15) 



*(J2) - ^ 2 ) exp (^-^!))*(#). (A16) 

Note that in rotated coordinates R = (i?i,i?2) — ► -R' = 
(i?2j -fti)- The electric field transforms as a vector under 
reflections, and the magnetic flux as a pseudoscalar. 



Time reversal. Time-reversal symmetry is imple- 
mented by the antiunitary operation 



fra > £rf ra , 
fra * Crfra- 

The resulting continuum operation is 

[* t (<7*)(<M 3 )] T , 



(A17) 



(A18) 



[(ir 3 )(i/z 3 )*] 



The electric field is odd under time reversal, while the 
magnetic flux does not transform. This reverses the more 
familiar situation of real electromagnetism, where elec- 
tric charge is invariant under time reversal but magnetic 
charge is odd. 



APPENDIX B: LARGE- TV/ DIAGRAMMATICS 
AND RG 

The starting point for the large- Nf expansion is simply 
naive perturbation theory in the gauge interaction vertex. 
The fermion propagator is 



1 



A 

fc 2 ' 



where we have introduced the notation ft 
bare photon propagator takes the form 



m q v 



<W + ( 



1 



, 1n<b> 

<i 2 



(9 2 /e 2 ) 



(Bl) 



The 



(B2) 



The rather unusual momentum dependence of the nu- 
merator is due to our choice of a "nonlocal" gauge - 
this choice is made purely for technical convenience as 
it results in a simpler form for the photon propagator at 
leading order in X/Nf. Finally we have the vertex 




(B3) 



The next step in constructing the large- Nf perturba- 
tion theory is to calculate the leading order photon prop- 
agator. Recalling that e 2 ~ it is easy to see that 
the leading contribution represented by the geometric se- 
ries 



(B4) 




Upon summing the series and taking the limit of small q, 
the full photon propagator is 

16 



(5, v + {i-l) q -^-)+0{l/Nj) (B5) 
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The perturbation series for any desired correlator is 
then easily built out of the 1/Nf photon propagator 
Eq. (|B5Jl . the fermion propagator, and the vertex. For 
example, the fermion Green's function (^'(k)^'(k')) — 
(27r) 3 (5(fc — k')G( 2 \k) is represented as 



■+0(1/N]). (B6) 



As discussed in Sec. lIIBl we can implement a renormal- 
ization group using the large- Nf expansion. As an exam- 
ple, consider the fermion Green's function and suppose 
we have added a single perturbation to the Lagrangian, 
represented by the coupling g. It is most convenient to 
explicitly keep track only of the anomalous part of the 
scaling - that is, we implicitly subtract off all contribu- 
tions to the Callan-Symanzik equation that give rise to 
the engineering dimensions of fields and coupling con- 
stants. This is denoted by writing primed versions of 
the appropriate quantities; for example, the engineering 
dimension of W is unity, so, denoting the scaling dimen- 
sion of ^ by A<f , we write A3, = 1 + A^. The Callan- 
Symanzik equation then takes the form 



+ 2A i + (t)'|] G(2, w = °- < B7 > 



APPENDIX C: IRRELEVANCE OF THE 
VELOCITY ANISOTROPY 

In this appendix we show that the velocity anisotropy 
for the fermions is irrelevant at the algebraic spin liq- 
uid fixed point, at least to leading order in 1/Nf. We 
do this by calculating the coefficient of 5/Nf in d8/M. 
In Refs. |52 and |U the same RG flow was calculated to 
leading order in 1/Nf but for arbitrary 5. Our calcula- 
tion is essentially equivalent to that of Refs. 112 andES 
the only difference is that we are interested here only in 
local stability and can work perturbatively in 8 from the 
beginning, which simplifies some technical aspects and 
makes the presentation more striaightforward. Our re- 
sults are in complete agreement with those of Refs. I22I 
and[13 

We employ the renormalization group approach dis- 
cussed in Sec. Ill Bl We need to calculate the Green's 
function G^ 2 \k), keeping terms of order <5, 1/Nf, and 
5/Nf, and then apply the appropriate Callan-Symanzik 
equation to determine db/dl. For ease of presentation we 
work in Feynman gauge (£ = 1). We have also carried out 
these calculations in an arbitrary covariant gauge with no 
effect on the final result. 

The anisotropy was discussed in Sec. 



IrnBI and the 

perturbation to the Lagrangian can be written 



K a = -iS^fi^^df, + ia^V, 



(CI) 



where we have introduced the notation 7^ = 71^.1 — 
72^,2- This term is represented by two vertices. The 



first is the correction to the fermion kinetic energy: 

» H » =-S^{k 111 ~k 2l2 )^-S^^. (C2) 

Here (i = (^171 — £272)- The second is the correction to 
the vertex: 




(C3) 



We need to calculate the order 1/Nf contributions to 
the fermion self-energy, up through linear order in 5. 
These are given by a sum of four diagrams 



S( 1 / JV /)(fc) = ^S i (fc). 



Here Sq is the isotropic (i.e. S — 0) contribution 



(C4) 



£o(fc) 



1G 



k+q 



d 3 q 1 r 1 



N f J (2tt)3 \q\ l^(fi+ 



:7m 



(C5) 

The other three diagrams are the anisotropy contribu- 
tions. We suppress the momentum labels since the mo- 
mentum structure is the same as Eq. (|C5|) . We have 



Si(fc) 



s 2 (fc) 



s 3 (fc) 



168fi 3 f d 3 q 1 r 1 



U&+ 4) 



(C6) 
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d 3 q 1 



7m 
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7m 



(#+ 4) 



(#+ 4) 



7m 



(C7) 



(C8) 



7m 



We evaluate these integrals using dimensional regular- 
ization, which introduces a mass scale /j, that roughly 
plays the role of a UV cutoff. Keeping track only of the 
logarithmically divergent parts, the results are 



3ir 2 N f 



*MI*I/m) 



(C9) 



and 

E 1 (fc) + E 2 (fc) + S 3 (fc) 



1 59 A// 3 - 

_^-^ln(|fe|/ M ). (CIO) 
15n z N f 



We can now calculate (k) to the appropriate or- 
der and use the following Callan-Symanzik equation to 
determine the flow of the anistropy: 



'd5\' d 
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FIG. 3: Representation of an arbitrary diagram contribut- 
ing to the two-point function of a fermion bilinear. The two 
dark circles represent the bilinear's matrix structure and are 
positioned at the points where the bilinear is inserted into 
the diagram. The shaded region is built from the ingredients 
of the large- Nf perturbation theory described in Sec. Ill Bl - 
fermion and photon lines, and the gauge vertex. The pairs 
(a,i) denote the SU(Nf) flavor (a = l,...,Nf) and Dirac 
spin (i — 1,2) indices of the fermions. 



The diagrams contributing to the Green's function are 



► a ► 



+ — ► » ► & > i ► 





and therefore 



G(2) (fc) = 12 



1 



ln(|fc|/ M ) 



(C12) 



(C13) 
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f f6 I 15ir 2 N f 



ln(|fe|/ M ) 



Noting that (dS/dl) 1 = (d5/d£) because S is marginal 
at infinite Nf, it is straightforward to apply Eq. (|Cllll to 
obtain = -4/3ir 2 N f and 



G4 



d6 



di 5ir 2 N f 



5. 



(C14) 



APPENDIX D: TWO-POINT CORRELATIONS 
OF N a AND M TO ALL ORDERS IN 1/N 

In this appendix we show that the fermion bilinears 
N a = -i^>T a ^ and M = —WW have the same scaling 



dimension to all orders in the 1 /Nf expansion. It is sim- 
plest to do this by a direct consideration of the two-point 
functions of these operators, defined by 



C M (x) = ((§tf)(a:)(**)(0)), 
C#(x) = ((*T a ^)(a;)(*T h *)(0)). 

Precisely, we shall show that 

Ctf(x) = X6 ab C M (x), 



(Dl) 
(D2) 



(D3) 



where A is an unimportant proportionality constant. 

Diagrammatically, all contributions to both these cor- 
relators can be represented in the form shown in Fig. [3] 
In the case of C$ the left and right dark circles repre- 
sent the flavor matrices T a and T b , respectively, while for 
Cm they are simply the identity matrix. The important 
point is that in both cases these matrices are trivial in 
the Dirac 7-matrix space. 

The shaded region represents an arbitrary combina- 
tion of the elements of the large- Nf perturbation theory. 
These are all trivial in flavor indices, so by SU (Nf) sym- 
metry the shaded region can only give a contribution pro- 
portional to either <5 Q/ 3<5 Q '/3< or 8 aa '5i3i3'- Let us consider 
these two possibilities in turn. 

The first of these corresponds to the a fermion-line 
extending through the shaded area and eventually joining 
onto the (3 line on the same side, and similarly for the a' 
and P' lines. So each of the external dark circles lives on 
a separate closed fermion loop, which is decorated with 
photon lines inside the shaded area. However, because 
the fermion bilinears are trivial in the Dirac space, it is 
easy to see that these fermion loops necessarily involve a 
trace over an odd number of 7 matrices, which vanishes. 
Therefore there are no contributions of this form. 

Next we consider the second possibility, where the a 
fermion-line extends into the shaded area and emerges 
on the other side to join the a' line, and similarly for the 
[3 and /?' fermion lines. In this case the two dark circles 
reside on the same closed fermion loop, which will involve 
an even number of 7 matrices and can be nonzero. Every 
diagram contributing to either correlation function will 
be of this form. Furthermore, for every diagram con- 
tributing to Cm, there is a unique diagram contributing 
to Off that differs only by the structure of the trace over 
this one-fermion loop. In the case of Cm this trace will 
have the form Tr(FiF2), where Ti is some matrix trivial 
in the flavor space. For Cff we have instead 



Tr(T a T 1 T b T 2 ) = Tr(T a T b )Tr{T 1 T 2 ) = X6 ab Tr^r 2 ), 

(D4) 

where A is a constant that is chosen once and for all by 
fixing a normalization for the SU(Nf) generators. This 
means that Eq. (|D3(I holds diagram-by-diagram, so N a 
and M have the same dimension to all orders in 1/Nf. 
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